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Abstract. We present the detection of pulsed gamma-ray emission from the Crab pulsar 
above 100 GeV with the VERITAS array of atmospheric Cherenkov telescopes [T] . Gamma- 
ray emission at theses energies was not expected in pulsar models. The detection of pulsed 
emission above 100 GeV and the absence of an exponential cutoff makes it unlikely that curvature 
radiation is the primary production mechanism of gamma rays at these energies. 



1. Introduction 

One of the most powerful pulsars in gamma rays is the Crab pulsar [21 [3], PSR J0534+220, 
which is the remnant of a historical supernova that was observed in 1054 A.D. It is located at a 
distance of 6500 light years, has a rotation period of ~33 ms, a spin-down power of 4.6 x 10 38 erg 
s and a surface mag netic field of 3.78 x 10 12 G @]. 

Within the corotating magnetosphere, charged particles are accelerated to relativistic energies 
and emit non-thermal radiation from radio waves through gamma rays. In general, gamma-ray 
pulsars exhibit a break in the spectrum between a few hundred MeV and a few GeV. Mapping the 
cut-off can help to constrain the geometry of the acceleration region, the gamma-ray radiation 
mechanisms and the attenuation of gamma-rays. 

Although past measurements of the Crab pulsar spectrum are consistent with a power law 
with exponential cut-off, flux measurements above 10 GeV are systematically above the best-fit 
model, suggesting that the spectrum is indeed harder than a power law with exponential cut- 
off [31 [5]. However, the statistical uncertainty of the previous data was insufficient to allow a 
definite conclusion about the spectral shape. In this paper we summarise the recent detection 
of the Crab pulsar above 100 GeV with VERITAS. 



2. Observation and analysis 

VERITAS, the Very Energetic Radiation Imaging Telescope Array System, is an array of four 
12 m diameter imaging atmospheric Cherenkov telescopes located in southern Arizona, USA 
[6]. After evidence for pulsed emission was seen in 45 hours of data from the Crab pulsar that 
were recorded between 2007 and 2010, a deep 62-hour observation was carried out on the Crab 
pulsar between September 2010 and March 2011. The observations were made in wobble mode 
with an 0.5 degree offset. After eliminating data taken under variable or poor sky conditions or 
affected by technical problems, the total analysed data set comprises 107 hours of observations 
(97 hours dead-time corrected) carried out with all four telescopes. The data were taken with the 
standard VERITAS trigger setting, and analysed with the standard VERITAS analysis tools. 
Details about the analysis can be found in pQ. 
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Figure 1. VERITAS pulse profile of the Crab pulsar at > 120 GeV. The shaded histogram 
show the VERITAS data. The pulse profile is shown twice for clarity. The dashed horizontal 
line shows the background level estimated from data in the phase region between 0.43 and 0.94. 
The data above 100 MeV from the Fermi-LAT [3] are shown beneath the VERITAS profile. 
The vertical dashed lines in the panels mark the best-fit peak positions of PI and P2 in the 
VERITAS data. 




Figure 2. Enlarged view of the main pulse or Figure 3. Enlarged view of the interpulse or 
PI of the Crab pulsar pulse profile. See text. P2 of the Crab pulsar pulse profile. See text. 



Figure 4. Spectral energy distri- 
bution (SED) of the Crab pulsar in 
gamma rays. The bow tie and the 
enclosed dotted line give the statis- 
tical uncertainties and the best-fit 
power-law spectrum for the VERI- 
TAS data using a forward-folding 
method. The result of a fit of 
the VERITAS and Fermi-LAT data 
with a broken power law is given 
by the solid line and the result of a 
fit with a power-law spectrum mul- 
tiplied with an exponential cut-off 
is given by the dashed line. Below 
the SED we plot \ 2 values to vi- 
sualise the deviations of the best- 
fit parametrisation from the Fermi- 
LAT and VERITAS flux measure- 
ments. 
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3. Results 

3.1. Pulse Profile 

The phase-folded event distribution, hereafter pulse profile, of the selected VERITAS events is 
shown in Figure [TJ The most significant structures are two pulses with peak amplitudes at phase 
0.0 and phase 0.4. These coincide with the locations of the main pulse and interpulse, hereafter 
PI and P2, which are the two main features in the pulse profile of the Crab pulsar throughout 
the electromagnetic spectrum. In order to assess the significance of the pulsed emission, we 
use the H-Test [7]. The test result is 50, which translates into a statistical significance of 6.0 
standard deviations that pulsed emission is present in the data. 

The pulse profile has been characterised by an unbinned maximum-likelihood fit; see the solid 
black line in Figures [2] and In the fit, the pulses are modeled with Gaussian functions, and 
the background is determined from the events that fall between phases 0.5 and 0.9 in the pulse 
profile (referred to as the off-pulse region). 

The positions of PI and P2 in the VERITAS data are thus determined to lie at the phase 
values —0.0026 ± 0.0028 and 0.3978 ± 0.0020, respectively, and are shown by the vertical lines 
Figure 1. The full widths at half maximum (FWHM) of the fitted pulses are 0.0122 ±0.0035 and 
0.0267 ± 0.0052, respectively. The pulses are narrower than those measured by Fermi-LAT at 
100 MeV by a factor of two to three. The energy-dependent narrowing of the pulses is a strong 
probe of the nature of the magnetospheric particle acceleration region and can be used to shed 
some light on its geometry, electric field, and gamma-ray emission properties. 

3.2. Spectrum 

The gamma-ray spectrum above 100 GeV was measured by combining the signal regions around 
PI (phase -0.013 to 0.009) and P2 (phase 0.375 to 0.421). This can be considered a good 
approximation of the phase-averaged spectrum. However, the existence of a flux component 
that originates in the magnetosphere and is uniformly distributed in phase cannot be excluded 
and would be indistinguishable from the gamma-ray flux from the nebula. Figure H] shows 
the VERITAS phase-averaged spectrum together with measurements made with Fermi-LAT 
and MAGIC. In the energy range between 100 GeV and 400 GeV measured by VERITAS, 



the energy spectrum is well described by a power law dN/dE = A x (E/150 GeV) Q , with 
A = (4.2±0.6 stat +2.4 syst -1.4 syst ) x 1CT 11 TeV" 1 cm- 2 s^ 1 and a = -3.8±0.5 stat ±0.2 syst . 
The detection of pulsed gamma-ray emission between 200 GeV and 400 GeV, the highest energy 
flux point, is only possible if the emission region is at least 10 stellar radii from the star's surface 

Combining the VERITAS data with the Fermi-LAT data we can place a stringent constraint 
on the shape of the spectral turnover. The previously favoured spectral shape of the Crab pulsar 
above 1 GeV was an exponential cut-off dN/dE = A x (E/Eo) a exp(—E/Ec), which is a good 
parametrisation of the Fermi-LAT [3J and MAGIC [S] data. We note that the Fermi-LAT and 
MAGIC data can be equally well parametrised by a broken power law but those data are not 
sufficient to distinguish significantly between a broken power law and an exponential cut-off. 
The VERITAS data, on the other hand, clearly favour a broken power law as a parametrisation 
of the spectral shape. The fit of the VERITAS and Fermi-LAT data with a broken power law of 
the form A x (E/Eo) a /[l + (E/Eq) results in a x 2 value of 13.5 for 15 degrees of freedom 
with the fit parameters A = (1.45 ± 0.15 stat ) x 10" 5 TeV" 1 cm" 2 s -1 , E = 4.0 ± 0.5 stat GeV, 
a = — 1.96±0.02 s t a t and f3 = — 3.52±0.04 s ^ a ^ (see solid black line in Figure 2). A corresponding 
fit with a power law and an exponential cut-off yields a x 2 value of 66.8 for 16 degrees of freedom. 
The fit probability of 3.6 x 10 -8 derived from the x 2 value excludes the exponential cut-off as a 
viable parametrisation of the Crab pulsar spectrum. The conclusions do not change if systematic 
uncertainties in the energy scale of both experiments are taken into account. 

4. Discussion 

The detection of pulsed gamma-ray emission above 100 GeV provides strong constraints on the 
gamma-ray radiation mechanisms and the location of the acceleration regions. For example, 
the shape of the spectrum above the break cannot be attributed to curvature radiation because 
that would require an exponentially shaped cut-off. In addition, assuming a balance between 
acceleration gains and radiative losses by curvature radiation, the break in the gamma-ray 
spectrum is expected to be at E^ r = 24 GeV ?7 3//4 \/?! where r\ is the acceleration efficiency 
(j] < 1) and £ is the radius of curvature in units of the light-cylinder radius [5]. Though £ can 
be larger than one, only with an extremely large radius of curvature would it be possible to 
produce gamma-ray emission above 100 GeV with curvature radiation. It is, therefore, unlikely 
that curvature radiation is the dominant production mechanism of the observed gamma-ray 
emission above 100 GeV. Two possible interpretations are that either the entire gamma-ray 
production is dominated by one emission mechanism different from curvature radiation or that 
a second mechanism becomes dominant above the spectral break energy. 
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